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Exploiting the correspondence between Wigner distribution function and a frequency-resolved 
optical gating (FROG) measurement, we experimentally demonstrate existence of the chessboard¬ 
like interference patterns with a time-bandwidth product smaller than that of a transform-limited 
pulse in the phase space representation of compass states. Using superpositions of four electric pulses 
as realization of compass states, we have shown via direct measurements that displacements leading 
to orthogonal states can be smaller than limits set by uncertainty relations. In the experiment 
we observe an exactly chronocyclic correspondence to the sub-Planck structure in the interference 
pattern appearing for superposition of two Schodinger-cat-like states in a position-momentum phase 
space. 
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It is the superposition principle that leads to inter¬ 
ference and diffraction phenomena, determines evolu¬ 
tion of wavepackets in classical and quantum systems. 
When applying to macroscopic objects, the interpreta¬ 
tion paradox of wave function in superposition arises 
from Schrodinger’s gedanken experiment on cat states [l|. 
By a Schrodinger-cat-like state in quantum optics one 
usually understands a coherent superposition of two co¬ 
herent states, say, |o) + | —o). Experimentally, such super¬ 
positions were created, e.g., in the atomic and molecular 
systems 0,i, superconducting circuits ii, and quan¬ 
tum optical setups [H, 0] • R has been noted by Zurek [8| 
that a superposition of two Schrodinger-cat-like states 
(four coherent states in total, |o) + | - o) + \ia) + | - zo), a 
so called compass state) in a Wigner phase-space descrip¬ 
tion [il gives rise to the interference structure changing 
rapidly on an area smaller than a Planck’s constant h. 
The result is counterintuitive because the Heisenberg’s 
uncertainty principle sets a limitation on the simultane¬ 
ous resolutions in two conjugate observables. 

It turns out that the sub-Planck structure determines 
scales important to the distinguishability of quantum 
states [H , thereby, potentially has an impact on an ultra¬ 
sensitive quantum metrology Bini and could affect ef¬ 
ficient storage of quantum information MM- In prin¬ 
ciple, these sub -Planck structures could be used to im¬ 
prove the sensitivity in weak-force detection MM and 
to help maintaining a high fideliW in the continuous- 
variable teleportation protocols m M- In practice, 
the classical wave optics analogues of sub-Planck struc¬ 
tures in time-frequency domain were observed experi¬ 
mentally so far only for superposition of two Gaussian 
pulses [l8|, [l9| . Obviously, a single Schrodinger-cat-like 
state offers a sensitivity to the perturbations only in 
one direction: perpendicularly to the line joining the co¬ 
herent states. To provide sensitivity in all directions, 
another pair of coherent states is needed. Theoretical 
proposals for generation of compass states include in¬ 
teraction in cavity-QED systems [TsI, [ 2 ^, evolution in 
a Kerr medium 1^. 1^ and fractional revivals of molec¬ 
ular wavepackets |23|. Nevertheless, there was no ex¬ 


perimental proof in the more complex and demanding 
case of superposition of four pulses. The main technical 
challenge in the preparation of states separated simul¬ 
taneously in two conjugate coordinates comes from the 
difficulty in keeping the coherence among them. Further¬ 
more, performing measurements of such tiny phase space 
structures is far from trivial. 

In this Letter, we report a direct measurement of a 
compass state superposition of four optical pulses and the 
corresponding interference pattern in the time-frequency 
domain. The mathematical equivalence between the elec¬ 
tric field of an ultrashort pulse and the quantum me¬ 
chanical wave function of the same shape, enables us 
to observe the phase space structures through the time- 
dependent spectrum of light. The experimental data ob¬ 
tained from frequency-resolved optical gating (FROG) 
measurements of light pulses reveals sub-Planck structure 
analogues corresponding directly to the compass states. 
In the interference patterns, areas smaller than that of 
transform-limited pulses are measured, illustrating that 
displacements leading to orthogonality of compass states 
can be smaller than the Fourier limit imposed on the 
pulses forming these superpositions. 

As illustrated in Fig.[T]^b), in our second-harmonic gen¬ 
eration (SHG) FROG measurement the input optical 
pulse E{t) is split into two replicas with a tunable time 
delay by passing through a standard Michelson interfer¬ 
ometer. These two mutually delayed replicas mix in a 
nonlinear crystal such as the Barium borate (BBO) 
crystal used in our experiment. Then, a spectrally re¬ 
solved sum frequency signal is recorded for each time de¬ 
lay r. The resulting time-frequency map (spectrogram) 
implemented by the nonlinear self-optical gating mecha¬ 
nism can be formulated as the following function for an 
initial field E{t): 
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For pulses with real envelopes and a linear phase, the 
spectrogram m is easily mapped to the Wigner quasi- 
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probability distribution W{q,p) [9[: 
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as 
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Here, F*{q) denotes the complex conjugate of F{q), h is 
set to 1, and our time and frequency domains are rep¬ 
resented by (r, cj). Moreover, under these conditions, all 
cross-sections of the FROG spectrogram correspond to a 
scalar product between the probe state E(t) and its ‘twin’ 
appropriately shifted in time or frequency [l8|,[^. Thus, 
zeros of the cross-sections appear for these values of time 
and frequency shifts that lead to orthogonal states. In 
the following, we check the scale of the phase space dis¬ 
placements resulting in distinguishable states based on 
this property. Let us stress that the data presented here 
comes from direct spectrogram measurements and, un¬ 
like typical FROG applications, does not involve steps of 
state reconstruction. 

In our experiment, an Er-doped fiber laser with 
1564 nm center wavelength, 37 nm spectral bandwidth 
and 5.68 MHz repetition rate is used as the light source. 
To construct a compass state, i.e., the superposition of 
four coherent states, a pulse shaper is introduced to split 
the input pulse in both time and frequency domains, 
as illustrated in FigUJa) [26|. Here, a telescope is em¬ 
ployed to improve the spectral resolution of the pulse 
shaper by expanding the beam diameter from 2.8 mm 
(at the collimator) to 4.5 mm. The input pulse is directly 
shaped in frequency domain, though a set of grating, lens, 
and a spatial light modulator (SLM). The whole spec¬ 
trum occupies 320 pixels in the SLM, and the through¬ 
put of our pulse shaper is around 30%. The spectral 
separation is realized by blocking the central range of 
the input spectrum, see Fig. 2(a); while the temporal 
separation 2to is achieved by imposing an extra mask 
function Mslm(<^) = cos{ujto)e~^^^^ via the SLM [^ . 
The spectrally and temporally separated coherent pulses 
are further made transform limited by compensating the 
residual spectral phase via the phase modulation func¬ 
tion of the pulse shaper. The power spectrum and tem¬ 
poral intensity of an example are shown in Figs. Ha) 
and Hb), respectively. In total, an initial laser spectra 
was divided into four pulses separated simultaneously by 
2co’o in frequency and 2to in time, i.e., shaped into our 
time-frequency representation of the compass state. For 
practical reasons, in the experiment separation between 
pulses in frequency was kept constant while time separa¬ 
tion 2to varied from 1.5 ps to 5ps. 

In the ideal case, when the pulses cut from the initial 
laser spectrum have Gaussian envelopes, the state gener¬ 
ated in a pulse shaper has the form 

Ein(t) ^ ^-(t-\-to)‘^-i(uj-uJo)t 



FIG. 1: (Color online) Schematic view of our experimental 
setup, including (a) Pulse shaper composited of a set of grat¬ 
ing, lens, and spatial light modulator (SLM); (b) SHG FROG 
composited by a beam splitter (BS) and a motorized stage in 
one arm to control time delay r, a nonlinear crystal (BBO) for 
the second harmonic generation (SHG), and a spectrometer. 
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FIG. 2: Typical experimental profiles for our compass states 
are shown as a function of (a) relative frequency Ur = uj/{27r) 
and (b) time t. Here, the four pulses are concurrently sep¬ 
arated by ojo/tt = 3.3 THz in frequency and 2to = 4ps in 
time. 


where 2to and 2uJo are the time and frequency separa¬ 
tions between the Gaussian peaks. For the clarity of 
notation, normalization coefficients and parameters de¬ 
noting width of the wavepackets are omitted in Eq. 
Depending on the context, in the text that follows we will 
use either a regular frequency u or the corresponding an¬ 
gular frequency uj = 27ny. The compass state from 
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FIG. 3: (Color online) Compass state in the phase space ob¬ 
tained by (a) FROG map measured for to = 2ps; and (b) 
numerical simulation calculated for a superposition of four 
perfect identical Gaussian pulses. 
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Eq. dD) is constructed as a superposition of two mutually 
delayed pairs of pulses with the same carrier frequency 
in each pair; it corresponds to the original compass state 
from Zurek paper [H rotated by 7r/4 in the phase space. 

A comparison between a theoretical simulation and ex¬ 
perimentally obtained time and frequency phase space 
maps of compass states is presented in Fig.[3l where spec¬ 
trograms corresponding to the same time and frequency 
separation between the cat states are plotted. Exper¬ 
imental data from SHG EROG spectrogram measured 
for to = 2ps is shown in Eig. [3fa); while the theoretical 
plot calculated analytically for perfect Gaussian pulses 
is presented in Eig. [3fb). In addition to four peaks rep¬ 
resenting Gaussian wavepackets and located at the four 
corners of the plots, characteristic patterns of interference 
fringes are clearly visible between every two of the peaks. 
Moreover, in the center of four perimeters, a chessboard¬ 
like interference pattern appears. In the following, we 
will verify that for large enough separation between the 
peaks this chessboard-like pattern is constructed from 
areas smaller than AtAuj = 1/2 indicated by uncer¬ 
tainty relation, which would correspond to sub-Planck 
areas in the position-momentum phase space. Eor the 
shape of our pulses is obtained by cutting the laser spec¬ 
tra, a slight discrepancy between the patterns formed by 
the real experimental profiles and theoretically consid¬ 
ered ideal Gaussian pulses can be seen between Eigs.[3l^a) 
and [3Kb). 

In the Wigner representation of a compass state, the 
middle interference pattern is build up from small rect¬ 
angles of alternating positive and negative values of the 


FIG. 4: (Golor online) Zooms of the central interference 
structure of FROG maps measured for (a) to = 1.25ps, (b) 
to = 1.75ps, and (c) to = 2.5ps. As a comparison, the corre¬ 
sponding simulations calculated for compass states build from 
perfect Gaussian pulses are shown in (d), (e), and (f), respec¬ 
tively. 


function. Even though recordings in the EROG spec¬ 
trograms take on only non-negative values, areas of the 
above-mentioned rectangles remain the same, i.e., the 
distances between subsequent zero lines do not change. 
Eor a large enough separation distances between pulses 
forming the superposition, these areas grow smaller than 
the area unit defined by the elementary uncertainty re¬ 
lation. In quantum mechanical version, the area unit is 
equal to h/2; while in a chronocyclic phase space, it is 
simply 1/2. It is worth noting that even for smaller sepa¬ 
ration distances “sub-Eourier” areas might appear in the 
EROG maps, namely, the areas smaller than those cor¬ 
responding to dispersion AtAuj calculated for any from 
the pulses forming the compass state superposition. To 
demonstrate how the interference structure appearing in 
the middle of EROG maps changes with change of initial 
parameters, in Eig.jH zooms of the central parts of spec¬ 
trograms for different time separations between pair of 
pulses are presented. Eigure|4] shows the Zoom of the cen¬ 
tral interference patterns obtained from the EROG traces 
for (a) to = 1.25 ps, (b) to = 1.75 ps, and (c) to = 2.5 ps. 
The corresponding theoretical plots made for superposi¬ 
tions of four identical Gaussian pulses are shown in in 
Eigs. SKd), (e), and (f), respectively. The comparison 
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FIG. 5: (Color online) (a) A central frequency = 
383.36 THz) cross-section of the FROG map measured for 
to = 2.5 ps. Here, blue dots depict measured intensity values 
for a given time delay r introduced between two input state 
copies in the arms of FROG apparatus, (b) Average areas, 
AtAcj, of the rectangles appearing between the zero lines in 
the interference structure of the central part of FROG maps. 
Values depicted by red or gray dots correspond to two inde¬ 
pendent sets of experimental data, respectively. Region below 
the uncertainty relation limit of 0.5 is denoted in yellow. 

serves to illustrate that areas between zero lines of FROG 
maps indeed decrease with an increasing separation be¬ 
tween pair of pulses. Again, as due to imperfections, 
amplitudes of pulses used in the experiment were not the 
same and the resulting interference patterns measured 
are not symmetric in respect to the central wavelength 
line. 

Finally, let us examine the cross-section of the FROG 
map measured for to = 2.5 ps along the central wave¬ 
length A = 782 nm {u = 383.36 THz) presented in 
Fig. [5fa). As mentioned before, the cross-sections of 
FROG spectrograms give values of the scalar product 
of a probe field with its ideal copy, but shifted in phase 
space in a direction perpendicular to the direction defin¬ 


ing the cross-section. It is clearly seen that the values of 
the function plotted in Fig.[5fa) decrease to zero between 
the subsequent peaks. A time shift equal to one half of 
the distance between the zeros results in the superposi¬ 
tions orthogonal to the initial compass state. This result 
is complementary to the one reported in Ref. [HI , where 
zeros in cross-sections defined by r = 0 were shown. Fig¬ 
ure [H^b) demonstrates how average areas AtAuj of in¬ 
dividual rectangles forming the central interference pat¬ 
tern, change with the change of parameter to- The av¬ 
erage values of AtAuj plotted in Fig.[5fb) are calculated 
for two independently collected sets of data (depicted by 
gray or red dots, respectively). These two sets of data 
are measured for slightly different pulse shapes and differ¬ 
ently set time delay resolution, i.e., different values of the 
smallest stage-motor step. It is clearly seen that in both 
cases for larger separation distances, areas between zeros 
indeed reach below limit imposed by uncertainty relation. 
Exactly these areas determine a scale of smallest change 
ensuring distinguishability of the mutually shifted states. 
As a last point, we would like to underline that the uncer¬ 
tainty relation is not violated. What we prove here is that 
a sub-Fourier change of initial state (shift in frequency 
and/or time) is enough to produce state orthogonal to 
the initial superposition. In the quantum mechanics this 
leads to a perfect distinguishability of states. 

In summary, using correspondence between FROG 
maps in a chronocyclic phase space and the Wigner dis¬ 
tribution function, we have demonstrated existence of 
the interference structure changing on areas smaller than 
that of a transform-limited pulse. We have performed 
FROG measurements of compass states realized through 
superpositions of four light pulses constructed with a help 
of a spatial light modulator (SLM), by manipulating the 
phase difference in the spectrum. Interpretation of the 
FROG spectrograms as maps of the values of a scalar 
product between probe pulses and their shifted copies, 
allowed us to show explicitly that the scale of changes 
leading to orthogonality of states indeed goes below the 
Fourier limit if sufficient time and frequency separations 
are simultaneously kept between the input pulses. The 
sub-Planck structure of the interference pattern in phase 
space representation of compass states not only mani¬ 
fests itself as a generalized Schodinger’s cat state, but 
also provides the platform to exploit quantum metrology 
and quantum information processing through ultrafast 
optics. 
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